Testing the universality of the many-body metal-insulator transition by time 
evolution of a disordered one-dimensional ultracold fermionic gas 
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It is now possible to study experimentally the combined effect of disorder and interactions in cold 
atom physics. Motivated by these developments we investigate the dynamics around the metal- 
insulator transition (MIT) in a one-dimensional (ID) Fermi gas with short-range interactions in 
a quasiperiodic potential by the time-dependent density-matrix renormalization group (tDMRG) 
technique. By tuning disorder and interactions we study the MIT from the weakly to the strongly 
interacting limit. The MIT is not universal as time evolution, well described by a process of anoma- 
lous diffusion, depends qualitatively on the interaction strength. By using scaling ideas we relate the 
parameter that controls the diffusion process with the critical exponent that describes the divergence 
of the localization length. In the limit of strong interactions theoretical arguments suggest that the 
motion at the MIT tends to ballistic and critical exponents approach mean-field predictions. 
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Studies of the interplay of interactions and disorder 
have flourished in recent years jl-Sj . Reasons for this re- 
newed interest include cold atom experiments more 
quantitative numerical simulations [2^ and novel theoreti- 
cal techniques 0, 0] ■ Adding further appeal to this prob- 
lem, numerical results for interacting ID bosons in 
a disordered potential contradicts rigorous mathematical 
predictions Q. Reasons for these discrepancies are not 
yet well understood [5] . 

Here we address a related problem: the time evolution 
of ID fermions with short-range attractive interactions 
by tDMRG techniques. We choose tDMRG over other 
techniques because the range of sizes that can be accessed 
is much larger. This is key to minimize finite size effects 
that might obscure the occurrence of localization. We 
focus on dynamical properties as the time dependence of 
the atom distribution is a natural observable in cold atom 
experiments. Disorder is modelled by a quasiperiodic 
potential @ that can be implemented experimentally [l[ , 



V{n) — Acos(27ra;n 



(1) 



with uj irrational, 9 G [0, 27r), and A > 0. In the non- 
interacting limit a ID tight-binding model with this po- 
tential and a hopping parameter J = 1 undergoes a MIT 
at Ac = 2 As attractive interactions are turned on 
Ac decreases [1]. It is thus possible to study the role 
of interactions at the MIT from the weak to the strong 
coupling limit. We employ the term MIT instead of 
superconducting-insulator transition because according 
to [11 quasi-long-range order is already broken when the 
insulator transition occurs. 
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The model— We employ tDMRG ^ to study the dy- 
namics of the L-site spin- 1/2 Hubbard model, 



n = -j 



L-l 

E 



L-l 
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(2) 
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in which c^.o- annihilates an atom at site i in spin state 

= ^i.t + "-i.i U < is the on- 



site interaction, and V{i) is given by ([1} with uj = {Vb — 
l)/2. The angle is chosen so that V{i) is symmetric 
relative to the center of the system. 

The tDMRG provides an efficient way to simulate the 
time evolution of a wavefunction obtained with DMRG. 
Our initial configuration {t — 0) is the ground state of 
the Hamiltonian where the disordered potential ([Ij is 
replaced by a simple potential well of width £ = 64 and 
depth D = 10 centered at the origin, 

Vt<oii) = De{\x,\^£/2), (3) 

where Q is the Heaviside function and Xi = i — {L — l)/2 
is the location of the site relative to the center of the 
system. For < > we compute the real-time evolution 
(t > 0) of this ground state under the Hamiltonian TL 
for L = 256 after the potential well is replaced by the 
quasiperiodic potential H is broken into terms af- 
fecting only two neighbouring lattice sites. The time 
evolution operator e"*^'^*, decomposed using the sec- 
ond order Suzuki- Trotter breakup, is iteratively applied 
on the ground state obtained by finite system DMRG. 
The time step At, measured in units of h/J, satisfies 
0.01 < At < 0.05 and m — 200 states have been kept in 
the DMRG simulation unless noted otherwise. 

Before we proceed with the calculation we provide a 
brief overview of previous research on this model. In the 
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non-interacting limit, ^7 = 0, the MIT is described by 
a process of anomalous diffusion controlled by the 
multifractal dimensions of the spectrum The local- 
ization length ^ (X |A — Ac|~'' diverges at the transition 
with v 1 [llj. For A = the model is exactly solvable 
[l^ for all t/'s. For ^ 1 it is mapped onto a weakly 
interacting hard-core Bose gas with a rescaled hopping 
parameter J' w .P/\U\ [13]. This suggests that the MIT 
will occur at Ac ~ J'^/\U\. 

For finite disorder and interactions there are already 
several studies about the static properties of Eq. © 

[Tsj and related models (l6l - l20 | . The dynamics of an 
interacting ID Bose gas in a quasiperiodic potential was 
first investigated numerically in [16|. For a more recent 
study in which interactions are treated in a mean-field 
fashion we refer to [17]. In [1, [l^ it was found that, in 
a ID Fermi gas with attractive interactions, Ac = Ac(C/) 
depends on the interaction and that weak disorder can 
enhance superfluidity. Renormalization group techniques 
were employed in [20| to study the weak disorder limit of 
spinless fermions in the Fibonacci chain, a quasiperiodic 
potential that is critical for each value of the coupling 
constant. For sufficiently weak interactions it was found 
in [TJ] that the spectrum of the Fibonacci chain is still 
multifractal. However in the case of the potential ([T|) the 
system becomes an insulator for J7 < and A = 2. 
Results. — In order to investigate the dynamics of ^ we 
first compute the n-th order moment defined as 




U= -10 



X = 0.28 



U=-6 



^ = 0.25 



(4) 



in which is the many-body wavefunction at time 

t. Here, i = 0, — 1 runs over the site index and 



is the number operator at site 



We 

set L = 256 and the number of fermions per spin to 
N = 12. Initially fermions are confined to sites 96-159 
by the potential well Vt<o(i) ([3]). Then we study the 
time evolution after the potential well is removed and 
the quasipotential V{i) is switched on at t = 0. 

The results for and different A's are depicted 

in Fig. [TJ The values U — —6,-10 correspond to the 
regime of strong coupling where the interaction energy 
is larger than the kinetic and potential energy due to 
the quasiperiodic potential. We clearly observe in Fig. [1] 
arrest of diffusion for sufficiently large A. The critical dis- 
order Ac < 2 for which the MIT occurs decreases as |J7| 

increases. We have estimated Ac directly from {x'^{t)) by (i7o"rmal diffusion) 'g{L) a 7 
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FIG. 1. (Color online) U = -10: tDMRG calculation of 
{x^{t)) Eq. (gl) for, from top to bottom, A = 0.1,0.15,0.28. 
Diffusion is clearly suppressed at A — 0.28. The dashed line 
is the best, a ~ 1.54, fitting function ([5} around the MIT 
A = Ac « 0.15. U = -6: tDMRG calculation of {x^{t)) @ 
for, from top to bottom, A = 0.15,0.25,0.40. The dashed 
line is the best, a « 1.48, fitting function ^ around the 
MIT A = Ac « 0.25. In both figures L = 256 and = 
12. The maximum time tmax that we can explore is dictated 
either by the stability of the tDMRG simulation {U, A) = 
(—10,0.28), (—6,0.40) or by the growing importance of finite 
size effects for t > tmax in the rest of cases. 



where xq, to and a are fitting parameters. We note that 
this fitting function is only an educated guess. We choose 
it because, despite its simplicity, it led to a good descrip- 
tion of the data. Other functions recently used in the lit- 
erature [H were also tried but the fitting was qualitative 
worse. Results of the best fit (see Fig.[T]) are presented in 
Fig. [2] for different values of [/ at A ~ Ac. It is observed 
that a depends on U and it is different from the one for 
U = 0, a ~ 1 ~ 2dii where dn is the Hausdorff dimension 
of the spectrum [3] . Therefore strong interactions modify 
substantially the dynamics at the MIT. 

This is an important result. According to the one pa- 
rameter scaling theory ^IJ the parameter a is related to 
the critical exponent v that labels the universality class of 
the MIT. Therefore different a{U) at the MIT correspond 
to different universality classes. An important concept of 
this theory is the dimensionless conductance g = Et/S 
where Et, the Thouless energy, is the energy related to 
the typical time for a particle to cross a sample of size L 
and S is the mean level spacing. For a disordered metal 
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identifying a narrow region of A's for which the dynamics 
becomes substantially slower than in the metallic region 
and also by an explicit calculation of the participation 
ratio [8]. In the latter the critical Ac at which MIT oc- 
curs, for a fixed (L,C/), is identified as a maximum of 
the participation ratio as a function of A. We have also 
found that Ac docs not strongly depend on the filling fac- 
tor provided that the chemical potential is far from the 
band edge. 

In order to fit the numerical data we employ the ansatz. 



{x'it))^xiii + it/tor) 



(5) 



oo for L — ^ oo since 
Et oc 1/L^ and (S oc 1/ . Analogously for an insulator 
g{L) oc e~^/^ decays exponentially. A MIT is character- 
ized by a scale independent g{L) = g^. Two mechanisms 
can lead to this scale invariance: localization effects that 
slow down the motion (a;^(t)) cx at the MIT and a 
multifractal spectrum (7| , with the Hausdorff dimension 
dn, that induces an anomalous scaling of 5 oc 1/L^/'^". 
Based on these arguments it was predicted in |2l| that in 
Z? = 1 a MIT will occur provided that 2d}i = a. In the 
non-interacting limit this relation was verified in 0, [l3] • 
In the presence of repulsive interactions it has been 




FIG. 2. (Color online) a (see text around ([5]) for details) 
versus the critical exponent v. The latter is obtained 
by fitting to AS = ae"^^^-^"^" with a,h,v fitting 
parameters, L = f3 and lj = 5/f3 in ([T]). Fitting is 
restricted to a small region A > Ac such that ^ < L. 
Error bars were obtained by considering the stability of 
the results under small changes in Ac and the fitting inter- 
val. From left to right, the points correspond to (17, Ac) = 
(-f 0, 0.15), (-6, 0.25), (-3, 0.47), (-2, 0.70), (-1, 1.01), (0, 2). 
Qualitative agreement with the expression i' ~ (solid line) 
is observed. Inset: v and a (left axis) and 1/Ac (right axis) 
plotted against U . The dashed line is a guide to the eye. 



suggested [2^ the scaling theory must include two pa- 
rameters, g and the ratio between an energy related to 
interactions and the mean level spacing. For attrac- 
tive interactions, especially in a quasiperiodic potential, 
the situation is less clear. In any case the above argu- 
ments, together with the numerical results above, provide 
a rather compelling albeit qualitative picture of the role 
of interactions: as |J7| increases, a > 1 increases and the 
motion becomes superdiffusive. According to the scaling 
theory, the MIT can occur only if dn also increases. Phys- 
ically that means that interactions smooth out the fractal 
properties of the spectrum at the MIT. The smoothing 
will be substantial when the interacting energy is much 
larger than the typical size of the subbands induced by 
the quasiperiodic potential around the Fermi energy. In 
this large \U\ hmit, corresponding to hard-core bosons, 
the spectrum is no longer fractal ( c?h ~ 1) and there- 
fore the dynamics at the MIT a — 2(iH ~ 2 approaches 
the ballistic limit. The numerical findings of |14| and the 
semi-analytical results of [20| for spinless fermions fully 
support this picture. We note that Q many features of 
the many-body MIT are similar to those of a single par- 
ticle in a Cayley tree [23|. For this model a = 2 and 
V ^ 1/2 around the MIT. It is thus tempting to specu- 
late that these results also applies to the Hamiltonian ^ 
in the limit \U\ oo. 

Before we turn to the next observable a few comments 
are in order: a) the fitting interval is long enough for dis- 
order and interactions to strongly influence the motion, 
b) the motion is slower as \U\ increases. The length of 
the fitting interval (see below) increases accordingly. As 
a result, for \U\ ^ 1 the value of a is more dependent 
on the interval. It is thus likely that additional transient 
terms are present in ([S]) . We stick to ([5]) because the ad- 
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FIG. 3. (Color online) Participation number P{t) ([6]) for the 
same parameters as Fig. [T] A non increasing P(t) for t oo 
is a signature of localization. 



dition of more terms without a clear physical motivation 
would lead to ambiguous results, c) the maximum time 
that we represent in the figures, and that it is used in the 
fittings, was chosen so that both the numerical error ac- 
cumulation [t < 1000) and finite size effects that obscure 
localization are negligible. For the latter the maximum 
time strongly depends on [/, A. This maximum time ^max 
for which finite size effects are not important is chosen by 
imposing that the occupation number of the last five sites 
remains less than 0.01 and no sharp increases occurs for 
smaller times. For instance, around A « Ac, imax ~ 650 
for U = -10 but only t„iax « 75 for t/ = -1, d) the 
moments of the distribution, such as those depicted in 
Fig. [1] can be easily measured in cold atom experiments. 
Therefore this model is an ideal candidate for experimen- 
tal tests of the MIT in strongly interacting ID cold Fermi 
gases. 

In order to obtain information of the time evolution of 
the full many-body wavefunction we have also computed 
the time-dependent participation number [24] , 



(6) 



which, up to normalization factors, gives an estimation of 
the number of sites which, at a given time, are occupied 
(see 25| for more information). In an insulator P{t) will 
be constant for sufficiently long times but in a metal it 
will always increase with time. Even a steady increase 
indicates that at least some parts of the wavepacket can 
escape localization. Therefore P{t) is an indicator of lo- 
calization of the full wavepacket. In Fig. [3] we plot P{t) 
for U — —10,-6 and different A's. The results are fully 
consistent with the previous calculation of moments. The 
transition is located around the same Ac and no increase 
in time is observed for A 3> Ac. 

We now turn to further substantiate the non- 
universality of the MIT by an explicit calculation of the 
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critical exponent v. For that purpose we study the sen- 
sitivity of the ground state to a change of boundary con- 
ditions [III, 

AE^Ep- Ea (7) 

where Ep and E^ stand for the ground state for peri- 
odic and anti-periodic boundary conditions respectively. 
As the MIT is approached from the insulator side AE oc 
e~^/'5 with ^ c>c |A — Ad^"^. We exploit this relation to 
find I', with L — 13. In Fig. [2] we present results for J^(f7) 
for different a (t/) at A = Ac. It is observed that as U in- 
creases v decreases from its non-interacting value i> ^ I. 
This is an additional indication that the MIT in many- 
body systems is not universal. However the expected 
approach to the mean-field limit = 1/2 for J7 3> 1 
seems to be slow. Theoretical arguments [111, [2^ sug- 
gest that the anomalous diffusion, through a{U), at the 
MIT is directly related to the critical exponent I'iU) that 
labels the universality class of the MIT. The simplest ex- 
pression consistent with ideas and techniques employed 
in the non-interacting limit [ill 113, [26] is v = ^ , which is 
in qualitative agreement (see Fig. [5]) with the numerical 
results. Finally we note that the calculation of z/, a is 
rather crude and subjected to substantial uncertainties 
in the fitting procedure. This is especially true for the 



U = —6,-10 for which the value of a is rather sensitive 
to both the fitting interval and the details of the fitting 
function ([S]) . For instance for smaller intervals and fitting 
functions including additional transient terms the values 
of a tend to be larger. 

In conclusion, we have carried out a tDMRG study of 
the MIT in an interacting ID Fermi gas in a quasiperi- 
odic potential. The main results of the paper are: a) the 
dynamics around the MIT is well described by a process 
of super-diffusion, b) the MIT is not universal ~ critical 
exponents depend on the interaction strength and slowly 
approach mean-field predictions for sufficiently strong in- 
teractions -, c) based on scaling arguments '2ll| we pro- 
pose that for strong interactions the dynamics tends to 
ballistic and the localization length ^ diverges at the MIT 
as ^ (X |A — Ac|~'' with i/ « 1/2, d) our results can be 
tested experimentally in cold atom settings. 
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